



UNITED STATES 

NAVAL POSTGRADUATE SCHOOL 










/o®# 



*""»iiiim»""'^ 



THE TRANSIENT LOAD SUPPORTING 
CAPACITY OF FLUID FILMS 

(FIRST REPORT) 

. . . By . . . 

ERNEST K. GATCOMBE 
Associate Professor of Mechanical Engineering 







A Report to the Bureau of Ordnance 

upon an investigation conducted under 

Project Order No. 17708 



20 DECEMBER 1950 



Technical Report No. 2 



U. S. Naval Academy Press — 13-18-50 — 30 



f f 









THE TRANSIENT LOAD SUPPORTING 
CAPACITY OF FLUID FILMS 



B-rK. Gat comb© (1) 



A study is presented In this 
paper of the thickness and the load 
supporting capacity of fluid films 
entrapped between rotating circular 
discs, under transient conditions. 

Part one, (see Fig. ij is devoted to 
a theoretical analysis of the thick- 
ness and the load supporting capacity 
of fluid films under steady state 
conditions. Part two presents the results of an experimental 
investigation of these thicknesses and loads under transient 
conditions. 

Car ful Instrumentation of this very complex non-steady- 

state problem, has yielded solutions which have not been 

obtainable by theoretical studies, except for cases involving 

assumptions which greatly alter the true results. It Is foxtnd 

that the approach motion of the surfaces bounding the film is 

respon sible for a much greater part of its load stipporting 

(l) Associate Professor of Navel and *aarin© EngineexJing, The 
United States Naval Postgraduate School, Annapolis, Hd, 

Prom 8 portion of a research project on the stud"; of the 
load supporting action of fluid films sponsored by the 
Bureau of Ordnance, Washington, D. C. 




Fis I. Discs Posessinc Rotary 
But No Approach Motion With 
Entrapped Fluid Wedge. 
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action than is the steady state rotational laotxon. 

STATEMENT OF PROBLEMS 

The K8ln problem is the investigation of th* thlcknes*: 
end the corresponding load supporting capacity of fluid films 
entrapped between circular rotating discs on parallel shafts 
which have their center-to-center distance changed so aa to 
cause the film to become alternately thin the/i thick. 

METHOD ^ Aq‘TACK s 

An es’act analytical solution to the steady state problem 
was securod (see Part one). The author had previously obtamrri 
an approximate solution (2), 

The transient state problem was solved by experimental 
means 5 It io much too comolicated to solve by analytical 
methods if on® retains such factors as the deformation of th-» 
surfaces of the discs, surface roughness , .etc . These factor, 
occur in their m>t-ural way in the experimental goluticn. 
DEi!!CRJFl'IO :^ Of fE£ EXPERIMENTAL SET-UP (3) s 

Tso pa?.''il@l steel shafts (see Fig, 2B) wes;'*' 

•.fid ®ount€-d with adjuatable taper rollor bearings in haia^-.r-i 

ohe €t fitiieh attached to the upper beam with a 8tror..si;~T-»c • 

with an adjusting ssrsw for controlling the initial gap 

between diace. The other shaft was mounted in adjust . .1*3 

roller b earings in hangars secured to the lower . b@e.m which 
(Sy" See Bibliography 
(3) See Fig. lE 
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jeould b© vibrated tranayersely by means of a Westinghouse Fatigue 
motor. The frequency and amplitude of vibration were controlled 
by an oscillator and an amplifier. Ho lost motion was 
permitted anywhere, Immedlat© detection of lost motion was 
available. Each shaft was driven by a separate air turbine 
wheel mounted on the shaft, A D.C. excited coll waa pieced 

i 

in the space provided on the lower shaft. A pick-up coll was 
arranged on the upper shaft. These colls did not revolve with 
the shafts. As the lower beam vibrated a signal occurred in 
the pick-up coll. This signal was fed directly to a galvanometer 
type oscillograph fitted with an automatic camera which gave 
permanent records of the signal generated. Oscillograms 
were taken for eases of the discs revolving with no oil on 
their surfaces . a certain height of wave form on th© oscillo- 
gram results, Tlien, another oscillogram was secured without 
changing a setting of any kind, but pairmitting oil to flow 
onto .:h© surface ';f the lower discs on the oncoming sid«„ 

Th© second oscillogram showed a much reduced amplitude of 
vibration 'f the boem. The thickness of the film and its 
load c'luiu then be determined by methods outlined 
Several t-ets of ^such records were taken for different ^ < 

of the amplifier. 

Lead wires were rolled between the discs a 
number of times to Insure no '’spring beck” of the thickm®- 
of the leads with the lower beam supported by jacks tc 



prevent bending. The thickness of these leads, detex'in'tnect by 
interferometer methods, geve the neutrel gap thickness h.^ 
Other leads were rolled between the discs and pounded hj th® 
vibrator until the height of the scope wave form reached 
that rhicL it had when the gap contained air. These lead 
thicknesses gave the amplitude of vibration of the lower 
she ft discs for various power inputs to the vibration ri*- • 
collii and the pick-up excitatlon-coll , 

The pick-up epparstiis showed a nearly linear reap one r 
for thesf' small amplitudes of vibration (0,001 approx.). 

The dyiii-mlc spring constant of the beam was secured. Thu., 
there resulted a means of obtaining a continuous picture oi 
the thfckneass of the film and its load supporting capacity 
mder these tranaient conditions. 

The sxlstine experimental set-up could not be used 
to obtain a steady state solution of the problem, since the 
pick-up aoraratus did not respond well to A.G. excltat ■ 

A relc'..i'-5i:;' largr noutrel gap setting was chosen, w lfd., 
according to existing theory, would yield almost no ■: 
st?-tu ic.-'i '' cu; porting component. .Since there is i.' rl ' ' 'h' 
state 1 ■' id for the particular chosen neutral gap set, . ■ - 
beam may be ';,onsid6red as vibrating about its neutri'’ 
position. 
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UOMENCLATURB: 

* = pole distance 

a * abscissa of the bl-polar coordinate system, pressure 
angle 

p- * ordlnai-e of the bl-polar coordinate system 
f » stream function 

d « distance from the origin to the center of the disc in 
the bl-polar solution 

0 * U >/~2rho 

h® as reciprocal of "stretch factor" 
ho * neutral gap thickness between discs 
p * pressure 

P * load- supporting capacity of the film per unit face width, 
pitch point 

r « ^1^2 

T ^ V 
1 2 

rj^ and V 2 = radii of eurwature of the discs 
-1 * 

8 = tan - - — ±-— 

V 2 rho 

T » torque 

and Ug peripheral velocities of rotating discs 

U * + U2 

2 

u, V, and w » coordinate velocity components 
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Xq = effective width of band of contact 

__2 = Laplscian operator 

^ ' dx2 dy2 8*2 






viscosity at atmospheric pressure and stated temperature 



)x * viscosity 



CO « angular velocity 



p - density 



e 





+ 




J>C = hf by definition 



INTRODUCTION: 

For some time the author has been aware of the fact that 
the transient state of affairs, in many hydrodynamical problems, 
accounts for the relatively large load supoortlng capacity of 
fluid films. iiJost of our practical enpri'neerlnK applications. 
Involving fluid films such as the lubrication of various 
bearings end gear teeth, come under this category. However, 
these particular non-steady-state problems are, at present, wholly 
unsolvable by analytical means, at least in their entirety. 

It has been only recently that we at the Naval Post- 
graduate School have acquired accurate and delicate enough 
equipment to tindertake the solution of such a complex problem. 

Part two of this report presents the results of an investi- 
gation of the load supporting capacity of two oils, namely, 

Mavy SlOO and 2190T. The results have been limited by the 
capacity of our equipment, but they show clearly that 



relatively thin oil film loads and thicknesses can and have 
been measured by these methods. Modifications are already 
under way which will permit the Investigation of heavier 
film loads. 

We now focus our attention on the analytlcsl solution 
of the steady state problem. The bl-polar coordinate system 
is conveniently emoloyed, since it provides boundaries of 
exactly the nature required in this investigation. 

PROCEDURE : 

The general procedure followed in the solution of 
this problem is to ia) write the equations of viscous 
fluid motion in the osrteslan coordinate system, (b) , 
assume that the fluid is incompressible and that constant- 
vlscosity-steady-state-condltlons exist, (c) accomplish the 
transformation to the bl-polar coordinate system, id) obtain 
the expressions which indicate the characteristics of the 
fluid film. 

EQUATIONS OF MOTION FOR STEADY STATE CONSTANT VISCOSITY 
INCOMPRESSIBLE FLUID CONDITIONS 

Fig, 2 shows an elemental cube of f^e fluid with its 
indicated positive stresses. Texts 
on hydrodynamics present the follow- 
ing equations for viscous fluid 
motion in which the viscosity p. is 
considered constant. 

I • ~5I T 

• \ <«• 

• “inr r 

* EbKiiiiNTABr Cvmc uf Oil With Iwoicatkd Politiv*. 




8 



They are 

r Jfir 



ft-.. 1 



9-P J. U. f-/>y 5 ^'^ 
^ X 






if * i- a ' 



% 

W 



U) 



If the body forces X, Y, and Z. are neglected, the fluid 
assiiTned Incompressible, and steady state exists, then 
Equations (l) reduce to 



_=r ^ V^t-^ 

tfTX 









3s. 

■ 



— 



iA) 



^3 . 

^ ^ YW 
3Z 

If the equations of (2) are differentiated with respect to 
jf ^ and z re s oe c t i vely , and added, then the following 



equation results'* 

Z, 

Yt 



— o 



We are now concerned with two dimensional flow r ly,80 
Equation (3) may Le expressed as 



»■ ^±1 ^ 



Vr> = ( 



f 



i)r = ^ 



) 
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Equation (3) la now transformed 
into an equation in bi-polar 
coordinates. As has been 
previously mentioned, these 
coordinates permit boundary 
conditions which are exceedingly 
suitable to this problem. Ho 
approximations, in regard to 
film thicknesses, need be made. 

The equation which is used to 
transform from Cartesian coordinates 
to bl-polar coordinates is 




Pig. 3 Circular Discs 
Associated With The 
Bi-Poler Coordinate 
System 



X -^c - etj (5) 

Fig. 3 presents the cc, p plane of the bl-poler coordinate 
system superimposed upon the x, y, plane of the Cartesian 
system. The origin for both systems is at 0, The equation 

I 

of a circle which represents the circular disc of this 
problem is as follows^ 

X -i- i t! ~ ^ h. <K. ( 6 ) 

The distance (d) from the origin to the center of the circle 
Is given by 



d = a coth 



cc. 



(7) 
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Smell a ■repi'esfsn a trie distance from tlie oripln to the focal 
point of the Ll-p«.*lar s;/_tem. The raaius (rj^) of the disc ia 
expressed by 






= a csch tt. 



(8) 



Tne f' ‘ f -il • onshlps exist* 



rtl.i 3 



^ ah 01 



ly = 



cosh «c “ cos p ’ coa t cos p 

>e( ' jOndifiw, V of an Irfcornprsssible fluid wo may writ© 



..( 9 ) 



du 

dx 



dv 

dy 



We mav »1 wvile 



M- - 






0 






(10) 



3 y 



^ X ^ 









itMSKKcMMNZSCl 



X 



(iij 



eh ex© \ .5 th-:: - f r. ti- fuf tion. 

lii the 1 "rtesifj-i , we may write 






'9x^ 






In the bl-polai system Squatlon (12) Is written ss 

w L* f ( 

T ^ ( ^cc». "a/ssa.) 

g 

where (h) la the reclorocal of the 'stiorch factor’, 



M 



sj 
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It la expressed symbolically by 

/ 






(14) 



or 



h » 



cosh d - coa p 



(15) 



a 



The differential equation satisfied by is 



S' 15 s>y 

oa-^S* 



= O 



An appropriate solution is (4) 



..(16) 






osO 



•f*^^ Mh-i 

Vt= I I 



( Vv+ 0( ^ - 'X 



^in H 0^ C<e>5 '>^1 ^ 



where the coefficients A^j , , Cj^ , are arbitrary 

constants, if we write hH -X , the corresponding expression 
for may be determined with the aid of the differential 
eque tion 



(4) See Bibliography 



a.'S= (c<.sh<- Ct,s/3)( +■ 






/ / ^7C ^ 

‘ C Sin^ 



c^/32- 



) 



3/5 

1 + c C,oS -}- (!lc* J /C 

It l8 fo\md t^at 

X — (/4ft + CojJj o( i- ( Co -tJ>o ®( ) 3//1 A 

-f ( A j i~ l^/ -h Cf ^/n k ~t J>/ K ) 

. <&o 

-/- 2 Cos\\& fAyiOak (yi-rl)<^ i~ ^>1 C€>^/l /J'K 

Q^»hUlV+0°( -f-JJ’vn Sinii (y^-l) ^ 



-) 






Slnc^ til© quantities.? end -/'If are conjugate function, 
equations of the following type are to he sotlsfledl 

_ ^'S)'S (c. 

2>“( 

Thus an expression for p in the ct, |3 system is 

^ C- ZJ>o^f^k<Sm3-hZlBa+^i)i/S~ C^sh^i :s^^- -^ ) 

ff . 

^ I [ C Ah j - (k- / j ) L n k 

t%=t ^ 

/ [cn-^0(C^+l^,*,)-(^-i)lC„-, '• ‘ 

Equation r2l) represents the general expression for t^- 
pressure p. If we consider problems in which th* -re*«*''» > 
is single valued^then all terms posaeswing d or p al-fjt 
will not be retained. The coefficient (Bq + Dx)^oi p may .• 



set equal to zero. 






4j 



i 
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Bo 



~ 0 (22) 

If , up are the velocity components in the positive cc, 
and p directions, w© have 



dY 



= slnp 



dX . jr 

a 



(23) 



and 



^ df dX- 7 
= d5“" 



‘P 



a 



slnh e. 



(24) 



We let the angular velocities of the cy?,inders be represented 
by t4> and write the following equations which represent the 
boundary conditions! 

For d = +di, f = ?i, and “ a'^cosech ,..,(25) 

For <t = ^ = "'*^1» ~ slnh (-cCi) - a*^cosech(-ai) 

(26) 

(27) 



By symmetry for ot = 0, ¥ = 0, and 






da2 



= 0 



Thus the expression for ^of tquf.tion (19> may be reduced. In 
the case of our problem, to 

cC( cosh (X cosp ) + si nrmc + slnh2ct cosp 28) 

where 

CL ^ C^sh J^OC^ 

Cos A Sin 



/• — €L Cosli 

t 3 ^ 

- CL cu 

C|^- 



4- Sfntiic^ 



(29) 
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If thia almplifled expression Equation ^28)^ is snbatituted 

into Equation (18)^ the following expression^ appropriate to 
this problem^ may be written? 



CL^ ^ ‘ 4 ^ C| Cos/?— CoSyS Co$iio( 



The derivative with respect to a yields 



A “ 4 Cf [e^shcc c^sfi «* Z c&tp Cask Zc(-f- C&sh Zo( j--- (s/J 

Inteeratlnff partially with respect to p following the statement 
of Equation i 20) gives 



Q.f>= C't'f/'C, 



r 

€psh^ Sm/B 



■^Cashz^ I 



Now a combi tlon of facts, t hat there is symmetry, that the 
fluid is assumed ln''apable of exerting tension stresses, that 
we are dealing with conversing and diverging spaces between 
the discs, require that the pressure shall pass through zero 
at the orisln where p = tt. 

Tviis yields G ^ 0* Thus our final expression for the pressure 
p, in the€|^ p system, is 

Co5^ oC ^ CojrA ^,#’5 ) 

^here 

a ^ 

^C•s^ 3 A 




— - 
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An •zpresslon for tho maxlnum presauro (Pmix*)* along tho 



z-Azlt or a * 0, la 

„ 31/^ C, 



(36) 



An axpresslon la now obtalnad for tho load-earrylng 
eapaelty (f) par nnlt faeo width of tho fluid wodge. Vo note froa 
Fig. 3. that an ezproaalon for this foreo la 



JP ^ Cos <p d6 

o ' 

It la naeaaaary to oxpraaa p of Equation 36 In terma of 



romeaborlng that cr aay ho oonalderod constant for any glren 



alza disc. Expressions for alnp and slnS6 are tho following: 



^ sz (37) 

di-f-rcoscp 

5 /^ ^ /3 ^ ^(3Q) 

(<J -t r Cos 

An expression for P la 




. ....(41) 



•ad 2 cosJ> Z 




We further observe that nearly all the load, carried by the 
fluid wedge. Is distributed over a very swall area. Let us 
arbitrarily define the width of the supporting band of contact 
Xq as that distance from the origin, measured along the x-ojcIl, 
to the point where the pressure drops off to one per cent of 
Its maximum value* An approximate expression for Xq is 

» 6.3 a (42) 

From, previous work by this author (5) it was found that 
an approximate solution to this problem gave the following 
set of equations } 

For the presstire £ 




• •eeeeaoeeeewe 



(43) 



where 







(44) 




U Is tne peripheral velocity of the disc. 




For equal siae discs r ® 

25 




(6) See Bibliography 
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L4n) 

ii8) 

^Yzrh^ ' ion(8o*S^'^.^ { 4 ^) 

, . *^ . i ’ 

•© shall now apply these equations In anelyzlng a 
given problem. Since the exact solution is obtained for the 
case of symmetrical flow.only^ we shall take discs of equal 
size revolving at the seme speed CU , The discs are 
3.0CC0 Inches In diameter. Their aneuler velocity is (-183) 
redy^sec. The minimum thickness of the wedge is h^ = 0.00002*BI 
inches; The viscosity of the oil is JJ = 5.5 x lO”® lb sec/sq.1n. 
We shall compute the following by the exact and approximate 
methods and compare the results: 

1. The maximum pressure within the wedge 

2. The load-carrying capacity P of the wedge per inch 
of face width 

3. The width Xq of the supporting band of contact 

SOLITflOM OF PROBLEM; - EXACT METHOD • 

Equation (35) is the expression for the maximum' pres Sbre 

c 

sr— 

OJ 

^ C0sA 





For the maximum pressure Pmax 




For the load-carrying capacity P 



For the width 'of the band of contact 



Cl 

a 



-(SO 



To find we need to use the following i 
d-j^ = a coth ct]^ ; ^ = a csch (£i 



or ^ = cosh (Xi 

d = 1.5000 + O.OOOOlCBfc = 1.50001BBI inches 
so cosh = 1.000006667 
oCi = 0.00365 red. 
and a =y^sink(ix ~ 1.5000 x 0.00365 
a - 0.00547 



Thus 



XS.SX/a 



/ 



/, ®e«» 3 CS" 






^ -t Z^Qd o -psi 



Equation (39) is the expression for the load-carrying 
capacity P 




ZM Si/s A 

£L 



iK^'^ 






= -.-0.00084 

Kg = -5.97 

P = -226 lb in. of fece width (55) 



Equation (42) is an expression for the width of the supporting 
band of contact 




m 




1 



le 



Xq = 6.3a = 6.3 x 0,00547 « 0.0344 Inches ...,.(56) 

SOLUTION OF PROBLEM; - APPROXIMATE METHOD 
The tnaxlmum pressure Pmax. given hy 

TV’ 

t,„ = -4^ G — 4— ^ Ze^ 7«e 7*^* 

IQ f Q J ^ 

The load-carrying capacity J? is indicated by 

/» « 

JP “ 225 lbs. per inch of face width. 

The width of the band of contact is 

X - 0.0 3^' 

COHCLUSIOH CONCBRNINQ STEADY STATE SOLUTION 

W« find that the approximate solution derived earlier, gives 
results which are close to the exact answers. However, at this 
point, our results are of little account in the solution of the 
usual non-steady-state problem except for a guide in the solution 
of steady state problems. If we wish to get practical results 
for these problems in hydrodyaamlcs we are forced to get the 
solution of these non-steady-state problems through experlEisai* 
amalye^. W» n 0 w pro«#«d to Part two of this report » 



PART II 



INTRODUCTION : 

The thickness and the load supporting capacity of oil 
films entrapped between circular rotating discs on parallel shafts 
were investigated here by experimental methods. (see Fig. IE). 

The fluid wedge was formed by dragging the oil between the 
discs which were rotated in opposite directions. The wedge was 
made alternately thick and thin by vibrating the beam which 
supported the lower shaft. The D.C. excitation coil set up 
lines of flux which had to pass through the neutral gap which 
separated the surfaces of the two sets of double discs. The 
pick-up coll delivered a siRnel to the ctalvanometer oscillograph 
as the lower shaft was vibrated. An automatic camera was used 
to photograph the oscillograms. Oscillograms were taken with 
and without fluid wddges between the discs, xhe presence of 
• oil caused the beajn supporting the lower shaft to vibrate 
through rel*' -Ively small amplitudes. Small lead wires of the 
order of 0.020 Inches were rolled between the two sets of discs 
with the lower beam supoorted by Jacks to prevent bending. 

The thickness of the neutral gap hQ was measured by obtaining 
the thlckne.ss of the flattened lead wire. The neutral gap 
thickness ho was adjusted until the amplitude of the vibrating 
beam, for a power setting of the amplifier of 60, equaled 
the value of the gap hQ. 

The constant D.C. excitation current for the pick-up 
system came from two storage batteries. The D C. and A.C. 
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DOUBLE DISCS 




verious amplifier settings. Then the values of the height 
of the wave form on the scope were found for various settings 
of the amplifier with the angular velocities of the discs at 
zero. These values are plotted in <Jraph No. 1. Lead wires 
were then rolled between the discs, left in position with 
CO = 0, and the vibrator pounded them \intil the height of 

the scope wave form equaled that which existed for the case 
of air between the discs. Leads were taken and their 
thickness measured for ell six power settings of the amplifier^ 
except no. 50 where the leads became too thin to handle. 

EXPERIMENTAL RESULTS: 

It is clear now that at this point we have Information 
at hand for determining the magnification factor (M.F.) of 
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our set-uD. Por e'^emole. we find from Graph No. 1, that for 
an amplifier setting of 60 (without oil and ca .^ = 0) the 
lower disc vihrated through an amplitude of hQ = 0.0C1P7 
inches, while the height of the scope wave form v as 1.80 
inches. This 3-ives a M.F. = 914. Like\»^lse. when the amplifier 




was set at 30 the leads showed that the lower shaft vibrated 
through 8(3 = 0.0C099 inches, while the scop© gave ag = 

1.44 inch's. The M.F. for this setting was 1450. 

It was noticed that the Sg depended not only upon the 
amplitude of vi.br at ion of the disc , but also upon the 
speed with which the discs rotated. It thus became necessary 
to plot the scope ampl Itude against the amplifier settings 
for various values of the angular velocity . These 
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values are shown In Graph No. 1. The magnification factor 
M.P. will always h© determined with reference to the graph 
of CO = 0. The increase In the amplitude of the scope wave 

due to the rotation of the discs Is nearly nullified for the 
cases with and without oil, since the angular velocity was 
fixed . 

We Indicate now how one may use Graph No, l to determine 
the load carried by a certain thickness oil film wedge. Prom 
Table No. 1 we find the following: that the height of the 

scope wave form, corresponding to = 2400 RPM without oil 
with the amplifier at 60, 1s 2.08 Inches, that the amplitude 
of the scope wave form, corresponding to = 2400 RPM with 
N3100 oil with the amplifier at 60, Is 1,35 Inches. We note 
that 0,28 Inches of the scope amplitude are due to rotational 
effects, Jw the amplitude of vibration of the lower disc, 
without oil, with the amplifier at 60 and f = 425 cps Is 
Just sufficient to close the neutral gap setting of h^ = 

0.00197 Inches, Thus, for the case of no oil, with co = 2400 

RTM, with f = 425 cps, the M.P, Is found to be 1,80/ 0.00197 
or 914. The magnification factor Is defined as the ratio of a / 
a^j , This value 1s calculated with reference to the curve tw 
= 0, Por the case of the use of N3100 oil with = 2400 RPM, 
with f * 426 cps we find from Table one that the Sg = 1,35 
Inches, We recall that 0,28 Inches of this amplitude Is due 
to rotational effects, thus the value corrected to the c*i » O 
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TABLE No. 1 Nqv, L9_5T) 





^ = 2400 


RPM 


Temp. 70° F. 1 


No Oil 


2I9( 


rr ' : j 


1 


N 


5100 




AMPLIFIER 

SETTINGS 


Os INS. 


jO^INS. 


Ofl INS 


M.F. 


t INS. 


Fq lb. 


Fs ,. lb 


[Qs INS. 


Qd INS. 


M.F. 


t INS. 


Fo lb. 


F... lb. 


Po lb. 


60 


2.08 


TiT 


0.00127 


13 60 


000070 


140 ' 




^1.3^ 


0.00081 


914 


0.00116 


162 


73.3 


102.6 


50 


2.04 


182 


0.001 1 1 


1420 


0.00086 


140 




1.25 


0.00073 


1280 


000124 


158 


1 


98.4 


40 


1.88 


1.80 


0.001 10 


1420 


0 00087 


140 




1.17 


0 00063 


1268 


0.00134 


156 




93.6 


30 


1.70 


1.50 


0.00087 


1350 


0.00 MO 


1 17 




0.97 


000061 


1 156 


000136 


134 




81.9 


20 


1.05 


088 


0.00054 


1 146 


0.00143 


69 




0.63 


0.00045 


976 


0.00152 ' 


74 




48.9 


10 


0.45 


0 39 


0.00043 


635 


0 00154 


21 




0.36 


0.00035 


622 


0.00162 


22 




20.4 




^ 




1 oo = 1750 RPM 1 


No Oil 


1 2190 T 1 


1 




No Oil 


AMPLIFIER 

SETTINGS 


Q| INS. 


Os INS. 


Qd INS. 


M.F 


t INS. 


Fo lb. 


F,.lb 


0, INS. 


Qd INS 


M.F. 


t INS. 


F, lb 


^,1. “> 


p, lb. 


60 


2.00 


TIF 


0 00126 


1 341 


0.00071 


136 


15 


114 


0 00072 


1255 


0.00125 


164 ' 


53 3 


102.6 


50 


1.90 


1.69 


0 00100 


1445 


0.00097 


142 




1 12 


0 00069 


1252 


0.00128 


164 




98.4 


40 


1.85 


1.65 


0.00099 


1450 


0 00098 


132 




1.05 


0 00064 


1250 


0 00133 


154 




93.6 


30 


173 


1 30 


0.00076 


1309 


0 00120 


123 




0 89 


0.0005 8 


1214 


0 00139 


1 33 




81.9 


20 


1.00 


0.75 


0.00057 


1060 


0.00140 


68 




0 59 


0 00046 


976 


0 00151 


73 




48 9 


10 


0.40 ! 


0.30 


0.00043 


620 


0.00154 


20 




0.29 


000039 


564 


0.00158 


20 




20.4 








= 875 RPM 1 


No Oil 1 


1 2190 T 1 


1 N 3100 1 


No Oil 


AMPLIFIER 

SETTINGS 


0, INS. 


QgINS. 


Q ^ INS. 


MF 


t INS 


Fo 'b 


F,,lb. 


Q, INS 


INS. 


M F. 


t INS. 


F, lb 


F,,ib. 


p, lb. 


60 


180 


174 


0.00129 


1336 


0.00068 


136 


7 3 


1.07 


lO 00079 


1250 


0.001 18 


164 


26 7 


102.6 


50 


1.78 


168 


0.001 1 8 


1405 


0.00079 


134 




104 


'0.00077 


1320 


0.00120 


154 




98 4 


40 


1.72 


1.65 


0.001 1 7 


1404 


0 00080 


133 




0.95 


0 00070 


1262 


0 00127 


150 




93.6 


30 


144 


1.40 


0.001 10 


1320 


0.00087 


105 




0.82 


0.00064 


’ 1250 


0 00133 


154 




81.9 


20 


0.80 


0 84 


0.00064 


1250 


0.00133 


65 




0 55 


0.00047 


1046 


0 00150 


73 




48.9 


1 0 


0 33 


0 30 


0 00043 


620 


0 00154 


20 




0.26 


0.00037 


560 


0 00160 


20 




20.4 



DYNAMIC SP CONSTANT OF BEAM ^ 250,000 lbs / IN 
VOLTAGE TO VIBRATION MOTOR = 8V (ACCORDING TO LOAD). 

CURRENT TO VIBRATION MOTOR = 200 M A. (ACCOPOING TO LOAD). 

VOLTAGE TO VIBRATION MOTOR * 118 V D C. 

CURRENT TO VIBRATION MOTOR = 0 9 A. D C 
VOLTAGE TO PICK-UP SYSTEM * 12V. DC. 

CURRENT TO PICK-UP SYSTEM * 0 25 A DC 

t*ho-a<i minimum oil wedge thickness. 

TOTAL WT. on beam » I2.|| lbs 

NATURAL FREQ OF BEAM WITH WT = 12.11 lbs * 450 cp$ 

Fq * LOAD CARRYING CAPACITY OF FILM IN Ibs/IN. OF FACE WIDTH. 

F, , *LOAD CARRYING CAPACITY OF FILM ( FROM STEADY STATE SOLUTION) 
Pq = AMPLITUDE OF MOTOR OUT PUT FORCE. 



No Oil f ■ 425 opt 
w - 0 RPM 






. AmpL at 60^ 

' 






1 




NoOil f-425cps 

«4i ■ 0 






Ampl. at 30 




1 


Showing Decoy 
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curve, is 1,07 Inches. This gives a M.F. = 1326. This 
scope amplitude of 1.07 Inches on the u/ = 0 curve corresponds 
to an amplifier setting of 24, and thus a disc amplitude of 
0.00081 Inches. The amplitude of vibration of the beam was 
thus reduced from h© - 0.00197 Inches without oil, to 

0.00081 Inches with oil. We know now that the minimum oil 
film wedge thickness was 0.00116 Inches. The load on the 
film will be computed now. 

From the mechanics of the problem we observe that the 
beam may be regarded as massless with an equivalent weight 
placed at the center of the span. There is also the weight 
of the lower shaft, bearings, hangars, motor drive coil, etc. 
which may be considered concentrated at the midsection of 
the beam. This total weight = 12,11 lbs. We see that four 
forces act on the beam in the displaced position, namely; 
the beam restoring force F|j , the force of internal damping 
Fd , the force of the oil film wedge Fq, and the force of 
the vibration motor Fjjj. 

The amplitude and the 
phase relationship of these 
forces are shown in Pig. 3S. 

The sum of the vertical com- 
ponents may be written ass 




FiO, 3F 



mtu yo + Pq cos # - Fq - kjo = 0 



(57) 
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The value of Fq, the amplitude of the motor output force, 
was found for different amplifier settings, by placing SR4 
strain gages on a specially prepared drive rod. These values 
are given In Table No, 2, It was found that the amplitude 
of the motor output force F* depended upon the amplifier 
setting for a given frequency. Thus the motor output force 
was evaluated without the use of oil between the discs. 



The phase angle # Is shown to be 

P ^ SJL 
‘an ♦ = ber ' 



^ - '“'Vo.® 



(58) 



^n ) 

where b Is the damping constant of the system and Is evaluated 

from the die-away wave form shown In the preceding photograph. 

It was found that b = 0.95 lb. sec.y/ In. The value of the 

critical damping constant b^j, = ^ 4mk was found to 

be bgjp = 177 lb. sec.y' In. The weight on the beam was 

12.11 lbs. and the dynamic spring constant was 250,000 lbs./ in. 

The undamped natural frequency of the system was “ 

2830 rad. / sec. It was found that ^ = 5® - 10^. 

The amplitude of load Pq carried by the oil film wedge 

may now be fovind from Eq. 57, The disc or beam amplitude of 

Vibrations a^j or Jq was found to be 0,00081". The mass m on 

2 

the beam was 0,0314 , The frequency of vlbfation 

in. 



of the beam was 425 cps or 2670 rad, j sec. The value of 
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was 102.61 lbs. Thus the oil wedsre, the minimum thickness of 

which was C.CC116 ", carried a load of 

Pq = - 250,000 X 0.00081 + 102.6 x 0.006 + 09514 x 

2 

X 0.00081 (59) 

or 

Fq = 8ip0 lbs. /half inch face width (60) 

Fq = 162.0 lbs. /inch of face width (61) 

If one used the results of Part I (So. 56) to oredict the 
width of the supportlnp; band of contact Xq , then the average 



pressure within the film will be 



^avT 



lo 

A 



81.0 



= 4710 osi.(62) 



0.5 X 0.0344 

0ns may count on the maximiim pressure being approximately 
seven times the average. This would give a maximum presvsure 
of 



max 



= 33,000 psi 



(63) 



It will iiow be of considerable interest to compute the load 
which this same minimum thickness oil film wedge will support 
on the basis of steady state conditions. t<q . 39 or Be. 48 
may be used to obtain this result. In 8q. 48 
G = 12 ^ U \/~2rho 



//a = 1.29 X lO'^reyns ( ^'*3100 r 70°F) . 

U = 378" / sec. (2400 PPM - 3" di am. disc) 
ho= 0.00197 
G = 3130 Dsl 



( 64 ) 









I 




1 





Thus 
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^ - I 



V" 2rh<5 



» 8130 

-g — X 0.054 = 73.3 Ihe.^ in. of fac© width or 

F ss 36.6 lbs./ half inch of face width. 

This is the load which the oil wedge could be expected to 
carry on the basis of steady state conditions. It is 
observed that the non- steady state solution indicates that 
this same wedge can carry 81 lbs. The non-steady state 
solution then shows that the wedge can carry 2.22 times as 
much load as the steady state solution indicates. From an 
Inspection of Fig. 4B it will be observed further that the 
ratio of the largest value of Fq to the smallest value of 
Fgg is = 6. This means that the load supporting 

action of the film, resulting from the approach motion of 
the metal bcmdlng surfaces, can be six times as large as 
that resulting from the rotational effects of the discs. 

I In other words the non- steady state solution gives values 
which are six times as large as those of the steady state. 

This is a factor which we may no longer omit in our pre- 
dictions of the load carrying capacity of films. It is 
no wonder that we have not been able to account for the 

large loads which bearings and gears actually carry. Our 
theory is based on steady state conditions. 

The preceding point is the prominent part of this whole 
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research project . The oil film wedee in bearings and gears 
supports loads because they operate under non- steady state 
conditions . 




Information on use of strain gages j 

The amplitude of the motor output force was determined 
through the use of strain gages. Two A-7 strain gages were 
used on the flat section of the drive rod. A third tempera- 
ture compensating gage was employed. The strain indication 
was received on a cathode ray oscillograph, through the 
employment of a Baldwin Southwark strain calibration and 
amplifier instrument. Two gages were used so that bending 



of the drive rod could be detected. It was determined that 
a strain of 120 micro inches per inch corresponded to an 
amplitude ordinate height of C.7 inches on the cathode ray 
oscillograph with the gain as used. The cross-sectional 
area of the steel drive rod was A = 0.0100 sq. ins. 

TABLE II 



Amplifier 


Amplitude 


Ampl Itude 




settings 


Gage A in Inches 


Gage B in inches 


60 


1.99 


2.00 


102.6 


60 


1.90 


1.91 


98.4 


40 


1.80 


1.81 


93.6 


30 


1.58 


1.59 


81.9 


20 


0.95 


0.95 


48.9 


10 


0.35 


0.40 


20.4 



Information on oils used: 

N3100 - navy symbol : 
specific gravity at 60 / 60°F = 0.903 
A.P.I, degrees = 25,2 

Viscosity seconds, SUV, at 100®F = 1318 
n It II It i3o°F = 492 

» » It n 210®F = 90.9 

Viscosity index (Dean and Davis) = 72,9 

2190T- turbine oil 

specific gravity at 60/ 60°F = 0.8844 



A.P.I. degrees 


= 28.5 








Viscosity seconds. 


SUV, 


at 


100°F = 


420 


n « 


n 


ft 


130®F = 


193 


n It 


II 


n 


210°F = 


5G 



(Engineering Experiment Station Test 1 BDT031) 
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Additional Remarks 

During the operation of the test equipment. It was 
found that an optional method of measuring the value of hQ 
gave a good check on those values already found through the 
use of leads. All these test results ere reproducible. 

In the test progra m to date the value selected for 
the nelitral gap setflng (h© = 0.00197“) was relatively 
large so that the beam would. In effect, vibrate about Its 
static neutral equilibrium position, rather than about some 
displaced position. In future work. If small values of 
are used. It will be necessary to correct the film thickness 
to account for this fact. The material In the steel discs 
was SAB 1045 steel hardened to 45 RC. The surface roughness 
was 20 millionths of an Inch. 

The Irregular frequency components of vibration 
present In tlse various wave forms shown In the photograph 
result from various causes such as non-concentrlclty of 
the discs, non-unlformlty of the rollers In the bearings, 
and non-uniform speed of the discs. These components do 
not alter appreciably the results of the test. 

Conclusions : 

A prominent cheractorlstlc of the load supporting 
capacity of fluid films has been uncovered In this research work. 
It Is definitely shown that thin fluid films strongly 
resist being squeezed out from between metal surfaces. 
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especially if the bounding surfaces are a poroacbing each 
other at high velocities. In other words, the non-steady- 
state character of the motion of the boundlnsr surfaces of 
fluid films is responsible for their ability to support 
relatively great loads. 

Our steady state theory never has been able to 
account satisfactorily for the relatively la'ge loads 
which our bearings and gears carry every day. This theory 
indicates unreasonably thin film thicknesses for these 
members when used with heavy loads at low soeeds; many 
times these thicknesses are much smaller than the surface 
roughness of the metal bounding surfaces. Neither can 
the theory account for the loads which are carried 
successfully in very high soeed gearing such as in certain 
turbine gear units. In such cases the period of mating 
is very short, and the whole lubrication ohenomenon is 
non-steady-state, without question. It is highly probable 
that the prominent point uncovered in this work partly 
explains why such gears, actually v/ork. 

Two external tyoe cylindrical discs, which formed 
the metal bounding surfaces of the fluid film under investi- 
gation, were chosen in this work in preference to actual 
bearings, because of their basic character in relation to 
gears as w^ll as to bearings. The results of this work 
may be used in studies of gears as ’rell as in internal 
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and external bearings. For one of the cases examined in this 
work, the non-steady-stete figures collected from actual 
experimental tests indicated that the film under inveati- 
gation carried a load six times as large as the steady s^ate 
solution would have indicated. 

No attempt v/ill be made at present to formulate a 
non-steady- state theory to explain this non-steady-state 
phenomenon. But it xvil 1 not be difficult to do this once 
we have more test data to draw from in forming semi-empirical 
theories. We should consider this work a simple beginning 
in the investigation of a most complex problem. 

It has been the opinion of this author for some 
time that theory alone is an insufficient tool for com- 
pletely explaining some of these complex lubrication phenom- 
ena. Thus we at the Nsval Postgraduate School have soent 
the last three years in acquiring and testing accurate 
Instrximents for carrying out this mission. It cannot be 
overemphasized that we have simulated actual operating 
conditions. We have worked with actual steel discs with 
real lubricants. These discs are deformable; they have 
surface roughness; they have certain adhesive and cohesive 
properties. The lubricants act in their natural fashion. 
Above all we have worked with non- steady- state conditions. 

This problem should be pursued vigorously along 
these combined experimental and theoretical lines until 
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enough date is gathered so that we may formulate the proper 
semi-empirical theories to give important p. . Id* nee to 
designers in the fields of bearing and gears as well as in 
other related fields. 
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